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Introduction 

This  report  Is  a  summary  of  activities  related  to 
Project  Sherwood  at  the  Magneto-Fluid  Dynamics  Division 
of  the  Courant  Institute  of  Mathematical  Sciences  of  New 
York  University  from  January  I96I  to  December  I961.   Sub- 
stantial investigations  are  summarized  under  the  following 
subject  headings: 

Charged  Particle  Trajectories 
Cusped  Geometries 
Golllslonless  Shock  Theory 
Wave  Propagation 
Other  Subjects. 

There  is  also  a  brief  summary  of  personnel  information  and 

a  bibliography  of  reports  and  papers  published  in  the  period 

covered.   For  a  complete  earlier  bibliography  see  Sherwood 

Progress  Reports  No.  3,  NYO-2544  (1959).  and  No.  4, 

NYO-939^  (1961). 
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Personnel 

The  average  number  of  equivalent  full-time  scientific 
staff  members  was  approximately  l6,  accounted  for  by  approxi- 
mately 20  Individuals  who  spent  an  appreciable  amount  of  time 
on  this  project.   Approximately  11  man-years  were  contributed 
by  senior  staff  members  (Ph.D.  level).   In  addition  there 
were  3  assistant  research  scientists,  2  graduate  students, 
1  active  consultant,  and  approximately  2  clerical  and  adminis- 
trative persons. 
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Charged  Particle  Trajectories 

Adlabatlc : 

The  adlabatlc  motion  of  the  guiding  center  of  a  single 
charged  particle  In  an  electromagnetic  field  has  previously 
been  analyzed  by  Kruskal  and  others  as  an  asymptotic  series 
expansion  of  the  true  trajectory,  cf.  Berkowltz  and  Gardner 
[Ref.  1,  PR  4],   A  simple  penetrating  treatment  of  the  prob- 
lem, using  perturbation  methods  of  classical  analysis,  was 
given  by  Gardner  [Ref.  3,  PR  4],  yielding  a  very  easy 
demonstration  of  the  Invarlance  of  the  magnetic  moment  to  all 
orders  In  the  expansion  parameter.   Tanlutl  [13]  has  given  an 
Invariant  treatment  of  the  canonical  theory  of  the  motion, 
exhibiting  explicitly  the  Hamlltonlan  In  terms  of  the 
coordinates  of  gyration  and  the  drift. 

Non-adlabatlc : 


The  underlying  assumption  which  Is  required  for  the 
validity  of  the  adlabatlc  theory  Is  that  In  one  Larmor  gyration 
the  particle  experiences  only  a  relatively  small  change  In  the 
electromagnetic  field.   This  Is  certainly  not  the  case  In  the 
neighborhood  of  a  point  where  the  magnetic  field  Is  zero. 
General  estimates  of  non-adlabatlc  effects  are  not  available, 
nor  are  they  to  be  expected,  although  In  some  fortuitous 
situations  some  special  results  have  been  obtained  (by  Hartweck 
and  Schliiter,  and  by  Chandrasekhar )  .   A  detailed  numerical 
study  has  been  made  of  particle  orbits  In  a  magnetic  field  of 
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the  cusped  geometry  type;  see  van  Norton  [17],  and  Grad  and 
van  Norton  [ 21] .   Two  features  which  have  emerged  from  the 
computations  are  that  there  Is  a  narrow  transition  zone 
between  very  adlabatlc  and  very  non-adlabatlc  orbits,  and  that 
even  the  slightly  non-adlabatlc  orbits  are  so  complex  that 
their  description  Is  amenable  to  statistical  methods.   Satis- 
factory containment  Is  possible  even  In  the  absence  of  an 
adlabatlc  Invariant.   The  results  do  not  seem  strongly 
dependent  on  the  particular  magnetic  field  configuration, 
and  may  therefore  be  applied  to  other  problems,  e.g.,  mirror 
machine  containment. 

A  feature  of  this  class  of  slightly  non-adlabatlc  orbits 
that  may  yield  to  mathematical  analysis  Is  the  following 
phenomenon,  observed  from  the  numerical  experiments:   The 
guiding  center  remains  on  a  line  of  force  and  the  magnetic 
moment  remains  almost  constant,  except  at  the  point  where  the 
guiding  center  Is  closest  to  the  zero  point  of  the  magnetic 
field;  here  the  magnetic  moment  jumps  to  another  slightly 
different  almost  constant  value.   This  slight  change  In 
magnetic  moment  in  turn  produces  a  change  of  the  exact  point 
at  which  the  particle  is  reflected  at  the  mirror  and  this  in 
turn  produces  a  large  difference  in  phase.   Thus  two  particles 
initially  In  the  same  state  except  for  slightly  different 
phases  will,  by  the  second  time  they  pass  near  the  zero  point 
of  the  field,  have  vastly  different  phases.   Thus  one  can 
consider  the  successive  increments  in  magnetic  moment  as 
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essentially  uncorrelated  with  the  original  phase  of  the 
particle.  This  randomization  of  phase  makes  possible  a 
stochastic  treatment  of  this  non-adlabatlc  effect. 

Based  on  some  profound  mathematical  discoveries  of 
Moser  [33a] , Gardner  has  recently  shown  that,  essentially.  If 
the  magnetic  field  Is  strong  enough  a  particle  not  In  the 
adlabatlc  loss  cone  Is  contained  for  all  time .   This  result 
Is  not  the  outcome  of  approximate  asymptotic  analysis,  but 
stems  from  a  rigorous  mathematical  theorem  which  reveals  the 
complicated  structure  of  the  manifold  of  charged  particle 
orbits. 


Berkowltz  and  Gardner  [Ref.  1,  PR  ^] 

Gardner  [Ref.  3,  PR  ^] 

Tanlutl  [13] 

van  Norton  [ 17] 

Grad  and  van  Norton  [21] 

Moser  [33a] 


Cusped  Geometries 

Stimulated  In  part  by  the  Increased  experimental  activity 
In  cusped  geometry,  there  has  been  some  effort  to  review  and 
Improve  the  available  theory. 

Grad  [12]  presented  a  descriptive  survey  of  the  current 
theoretical  picture  of  plasma  containment  In  cusped  magnetic 
configurations,  together  with  some  Indications  and  cautions 
about  the  relevance  of  the  theory  to  points  of  contact  with 
experiments . 

Considerable  effort  was  made  to  refine  the  containment 
theory  and  general  results  were  obtained  for  cusped  plasma 
systems  In  the  situation  with  arbitrary  p  where  plasma  and 
field  are  mixed  together.   A  numerical  study  of  single  particle 
orbits  In  cusped  magnetic  fields  (van  Norton  [17])  was  used  to 
test  the  conjecture  (Grad  [Ref.  22,  PR  ^])  that  the  particles 
can  be  classified  Into  three  types:   (l)  those  whose  motion  Is 
adlabatlc,  like  particles  contained  In  a  mirror  machine,  (2) 
those  whose  motion  Is  wildly  non-adlabatlc ,  like  those  In  an 
Idealized  P  =  1  sharp  Interface  cusped  plasma,  and  finally  (3) 
a  relatively  thin  transition  zone  of  slightly  non-adlabatlc 
particles.   Thus,  to  obtain  estimates  for  containment  times, 
one  combines  the  mirror-type  loss  rates  for  the  adlabatlc  par- 
ticles, and  cusp-type  loss  rates  for  the  non-adlabatlc  ones. 
The  resulting  formulas  are  given  by  Grad  [12],  together  with 
modifications  caused  by  the  presence  of  electric  fields.   Since 
there  Is  different  scaling  with  regard  to  both  temperature  and 
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density  In  the  two  types  of  loss  mechanism,  considerable  care 
must  be  exercised  to  Interpret  containment  results.   Moreover, 
the  loss  rates  can  be  very  sensitive  to  electrostatic  potentials 
which  may  well  be  altered  by  slight  changes  In  experimental 
arrangements . 

In  developing  a  containment  theory  for  cusped  geometries, 
as  well  as  for  many  other  problems  In  plasma  physics.  It  Is 
Important  to  have  a  theory  of  the  structure  of  the  steady 
transition  layer  separating  a  uniform  vacuum  magnetic  field 
from  a  colllslonless  field-free  plasma.   A  solution  for  the 
thinnest  transition  layer,  corresponding  to  any  non-singular 
velocity  distribution  prescribed  deep  In  the  plasma,  provided 
there  are  no  trapped  particles,  was  given  by  Grad  [Ref .  25,  PR  4] 
A  striking  feature  of  the  solution  Is  that  the  magnetic  field 
penetrates  only  a  finite  distance  into  the  plasma,  whereas 
the  plasma  decays  exponentially  into  the  magnetic  field  region. 
The  scale  of  length  Is  the  Larmor  radius.   This  solution  was 
obtained  ignoring  the  electric  fields  arising  from  charge 
separation  in  the  case  of  particles  of  equal  masses.   Numerical 
computations  are  underway  (Weltzner,  Paskievici,  Sestero)  to 
take  these  electric  fields  into  account,  as  well  as  to  study 
the  effect  of  trapped  particles  on  the  transition  layer. 

The  thinnest  boundary  layer  is  desirable  for  stability 
and  low  particle  losses  for  cusped  geometries,  and  it  is  best 
as  well  to  minimize  the  cyclotron  radiation.   Using  the  sheath 
analysis  just  mentioned,  the  minimum  radiation  has  been 
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determined  by  Burkhardt  [22] j 

¥  =  1  X  10"-^^  T^  x?''^    , 

o   o    ' 

o 
where   W   Is  in  watts/cm.  ,  T   Is  in  electron  volts,  and 

3 
n   Is  the  number  of  electrons  per  cm.  .   This  amount  of 

radiation  turns  out  to  be  tolerable  for  plasmas  of  thermo- 
nuclear Interest,  but  this  is  the  minimum  that  can  be  achieved 
and  not  necessarily  a  good  estimate. 

Some  mathematical  analysis  is  being  carried  out  to 
determine  details  of  cusp  shapes. 


Grad  [12] 

Grad  [Ref .  22,  PR  4] 
Grad  [  Ref.  25,  Pr  4] 
Grad  and  van  Norton  [21 ] 
van  Norton  [17] 
Burkhardt  [22] 


-  9  - 


Colllslonless  Shock  Theory 

Investigations  continue  on  the  Irreversible  heating  of 
plasmas  by  shock-like  phenomena.   The  situation  treated  con- 
cerns plasmas  in  which  the  classical  collision  mechanism  for 
explaining  shock  structure  is  not  applicable. 

Morawetz  [3]  has  now  published  a  description  of  the 
internal  structure  of  steady  colllslonless  shocks.   This 
theory  has  now  been  modified  [30]  to  take  into  account  finite 
electron  pressures,  but  the  problem  for  equal  ion  and  electron 
pressures  remains  open. 

The  equations  to  be  solved  are  the  colllslonless,  steady 
Boltzmann  equations  for  ions  and  electrons  coupled  with 
Maxwell's  equations  for  the  fields;  a  self -consistent  solution 
is  obtained  based  on  asymptotic  expansions  in  the  lon-to- 
electron  mass  ratio.   The  electrons  are  treated  adlabatlcally . 
These  results  show  an  Irreversible  transition  from  a  constant 
state  in  front  to  an  oscillating  state  behind.   The  case 
treated  supposes  the  flow  is  one-dimensional;  and  at  x  =  -oo, 
there  is  no  transverse  magnetic  field,  the  ion  pressure  is  low, 
the  Alfven-Mach  number  is  roughly  less  than  2,  and  the  electron 
pressure  is  zero.   In  this  analysis,  the  scale  of  length  is 
the  geometric  mean  Larmor  radius,  and  in  any  such  case  the 
Irreversible  heating  goes  mainly  to  the  ions. 

The  development  of  methods  for  solving  self-consistent 
field  problems  is  of  much  broader  interest  than  merely  for 
shock  problems. 
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An  Important  advance  has  been  made  In  the  analysis  of 
strong  shocks  (Morton  [l8]),  using  an  adlabatlc  two-fluid  model 
of  a  collision-free  plasma.   Numerical  computations  of  the 
transient  behavior  of  a  plasma  when  a  perfectly  conducting 
piston  Is  accelerated  Into  It,  In  many  cases,  shows  that  the 
compression  waves  steepen  and  break.   This  led  to  the  discovery 
of  the  existence  of  steady-state  generalized  discontinuous 
solutions.   Moreover,  numerical  analysis  of  the  time -dependent 
problem  revealed  rapid  convergence  to  these  discontinuous 
solutions.   Although  these  shocks  give  an  irreversible  heating, 
the  relative  heating  of  the  Ions  and  electrons  can  only  be 
determined  after  a  detailed  analysis  of  the  Internal  structure 
of  those  shocks  Is  made. 

In  the  theories  described  above  all  quantities  are  func- 
tions of  one  space  variable   x  and  time   t,   and  the  magnetic 
field  B   Is  supposed  to  point  always  in  the  z-dlrection. 
Recent  work  by  Saffman  and  Gardner  has  shown  that  there  is  a 
much  wider  class  of  solutions  when  this  last  condition  is 
removed.   In  order  to  assess  the  significance  of  this  discovery, 
we  have  begun  numerical  solutions  of  time-dependent  problems 
with  transverse  magnetic  fields.   There  are  some  preliminary 
results  in  the  two-fluid  model.   For  example,  in  the  cold  plasma 
case,  with  the  magnetic  field  only  in  the  z-dlrectlon,  the  com- 
pression wave  is  seen  as  a  steep  rise  followed  by  an  oscillating 
tall,  whereas  in  some  oblique  field  cases,  there  is  an  oscil- 
lating precursor  followed  by  the  steep  rise.   The  feasibility 
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of  similar  comp-utatlons  In  the  Boltzmann  equation  case  is 
under  Investigation. 

Gardner  and  Morlkawa  [ 1  o]  have  demonstrated  that  the 
width  of  a  small  amplitude  single  pressure  pulse,  propagating 
normal  to  the  magnetic  field.  Is  larger  by  a  factor  approxl- 
mately   (1  +  ^k"  ^)     than  that  In  the  zero  temperature 
case.   Here  p   Is  the  ratio  of  gas  pressure  to  magnetic 
pressure,  and  m  :  m_   Is  the  lon-to-electron  mass  ratio. 
This  factor  Is  significant  unless  p   Is  extremely  small,  so 
that,  at  least  for  weak  shocks,  the  zero  temperature  approxi- 
mation Is  valid  only  for  p   very  small,  say  P  <  0.001. 


Morawetz  [3] 

Morawetz  [30] 

Morton  [ l8] 

Gardner  and  Morlkawa  [ lo] 

Grad  and  Blank  [9] 
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Wave  Propagation 

Liboff  [  16]  has  completed  an  analysis  of  the  propagation 
of  long  wavelength  disturbances  In  a  plasma  in  a  steady  mag- 
netic field,  using  a  kinetic  theory  with  a  collision  term. 
These  results  generalize  previous  work  of  Bernstein,  Allls, 
and  others.   The  analysis  Involves  two  linearized  Boltzmann 
equations  with  modified  collision  terms  which  are  coupled 
through  the  differences  in  temperatures  and  in  velocities  of 
the  electron  and  ion  gases.   This  analysis  catalogues  the 
infinity  of  modes.  Identifying  those  found  in  earlier  more 
special  treatments  as  well  as  discovering  new  modes;  It  also 
reveals  interconnections  between  certain  macroscopic  and  micro- 
scopic modes.   Many  of  these  were  not  previously  known,  and 
they  only  appear  when  the  theory  is  as  broad  as  this  one.   To 
cite  one  example  from  the  multitude  of  results,  there  is  des- 
cribed the  intricate  interplay  between  the  fast  and  Alfven 
modes  of  magneto-fluid  dynamics  and  the  kinetic  plasma- 
magnetic  modes. 

There  is  also  underway  another  approach  using  a  sim- 
plified version  of  the  Fokker-Planck  equation,  instead  of  the 
Boltzmann  equation  with  collision  term,  which  is  appropriate 
for  other  parameter  ranges. 

Weltzner  [26]  has  completed  a  study  on  the  classical 
problem  of  the  small  oscillations  of  a  plasma  In  a  constant 
state  without  magnetic  field.   At  time   t  =  0  a  certain  dis- 
turbance Is  introduced  at  the  origin  x  =  0,   and  the  subsequent 
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wave  motion  Is  followed.   In  contrast  with  earlier  results, 
the  Initial  value  problem  is  solved,  giving  much  more  infor- 
mation than  can  be  obtained  from  plane  wave  solutions  and 
dispersion  relations.   Quite  arbitrary  initial  distribution 
functions  for  the  plasma  can  be  prescribed  --  there  are  no 
unnatural  analyticlty  requirements  as  in  other  treatments. 
Moreover,  the  solution  can  be  obtained  for  three  space  dimen- 
sions.  The  solution  has  been  evaluated  in  various  asymptotic 
limits  for  large  time.   For  stable  initial  distributions,  one 
finds  naturally  Landau  damping,  and,  for  example,  the  decay 
of  the  electric  field  can  be  explicitly  determined  as 

E(x,t)-'  £21ist.  (^^g^^  t  +  |!^)  +  ...). 

yv  P    ^ 

For  unstable  initial  distributions,  such  as  for  the  two-stream 
instability,  exponentially  growing  waves  are  found,  and  their 
growth  rates  calculated.   In  some  cases,  the  actual  growth 
rate  is  smaller  than  that  which  would  have  been  predicted  from 
the  dispersion  relation. 


Liboff  [  16] 
Weitzner  [26] 
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other  Subjects 

Grad  [20]  presented  a  survey  of  plasma  theory.   This 
Includes  a  greatly  simplified  discussion  of  the  justification 
of  the  use  of  the  Boltzmann  equation  and  Fokker-Planck  equa- 
tion in  a  fully  ionized  plasma.   In  particular,  it  is  shown 
that  radiation  coupling  (cooperative  fluctuation  effects)  play 
essentially  no  role  to  "dominant"  order.   Also  included  is  a 
simplified  presentation  of  the  theory  of  transport  coefficients 
and  a  general  survey  of  the  methods  available  for  solution  of 
problems. 

Another  survey  (Grad  [23])  was  written  on  the  question  of 
the  choice  of  practical  systems  of  equations  to  be  used.  This 
includes  a  simplification  of  the  theory  of  a  guiding-center 

■X- 

plasma;   in  particular,  a  description  is  given  which  makes  no 
use  of  "first  order  drifts".   There  are  many  cases  in  which  it 
is  a  priori  clear  that  a  guiding-center  particle  treatment  can 
give  no  improvement  over  a  simpler  fluid  analysis.   Also,  it 
is  brought  out  that  the  best  available  theory  is  incapable  of 
predicting  more  than  rather  short  time  plasma  containment 
(on  the  basis  of  equilibrium  theory,  not  stability).   Whether 
this  is  an  actual  difficulty  or  a  weakness  in  theory  is  not 
yet  clear. 


Essentially  an  elaboration  of  H.  Grad,  "A  Guiding  Center 
Fluid",  AEG  Report  TID-7503,  February  1956. 
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Tanluti  [15]  investigated  the  effect  of  finite  conductivity 
on  the  propagation  of  the  magneto-fluid  dynamic  slow  wave. 


Grad  [20] 
Grad  [23] 
Taniuti  [I5] 
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